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Data on the thermophysical properties of NZP ceramics with the NaZr,(P0 4 ) 3 -type structure are summarized. 
A system of thennal expansion regularities is proposed making it possible to predict NZP materials with con¬ 
trollable thermal expansion, including ultralow expansion with near-zero anisotropy. The thermodynamic 
functions of the reactions of synthesis of NZP phosphates are calculated and the applicability of ceramic tech¬ 
nology is justified. It is proposed to use NZP ceramics as heat-insulating materials capable of working under 
abrupt temperature shifts. 


The NZP crystallographic family with the basic structure 
of phosphate NaZr 2 (P0 4 ) 3 includes compounds and solid so¬ 
lutions described by the chrystallochemical formula 

(M1) 0 ^i(M 2) 0 ^ 3 {[L 2 (XO 4 ) 3 F } 3oo , 

where {[L 2 (X0 4 ) 3 Y~ } 3 oo is the skeletal framework (p is the 
charge of the skeleton); (Ml) 0 ^ 3 , (M2) 0 ^. 3 are the types of 
non-skeletal cation positions denoting the number of posi¬ 
tions in each type. 

The skeletal framework is formed by L0 6 octahedra and 
X0 4 tetrahedra joined by their vertices (Fig. 1); the vacancies 
Ml and M2 have different sizes and are convenient for incor¬ 
porating a wide range of charge-compensating cations. 

The majority of the representatives of the NZP family 
contain phosphor as the anion-forming element X. There are 
NZP-structure phases known in which phosphor is fully or 
partially replaced by silicon, germanium, arsenic, vanadium, 
sulfur, molybdenum, or tungsten. The structure-forming cat¬ 
ions L with predominantly covalent metal - oxygen bonds 
have degrees of oxidation ranging from +1 to +5. These are, 
for instance, Nb 5+ , Ta 5+ , Ti 4+ , Hf 4 , Ge 4+ , Sn 4+ , Mo 4+ , U 4+ , 
Sc 3+ , Y 3+ , Ln 3+ (lanthanides), V 3+ , Cr 3+ , Fe 3+ , Al 3+ , Ga 3+ , 
ln 3+ , Ti 3+ , Mg 2+ , Mn 2+ , Cu 2+ , Co 2+ , Ni 2+ , Zn 2+ , Na + , and 
others. The positions L of the skeleton are frequently occu¬ 
pied by a combination of these cations. The non-skeleton 
positions Ml and M2 may be filled by the same or different 
cations, mainly low-charge ones with degree of oxidation 
from +1 to +4, or they may stay vacant depending on the de¬ 
gree of oxidation of the ions in positions L and X and the 
electroneutrality of the compound in general. 

The resistance of an NZP structure to fluctuations in the 
chemical composition, temperature, and pressure is deter- 
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mined by relatively rigid structural fragments L 7 (X0 4 ) 3 with 
strong bonds inside the group and uniformly distributed va¬ 
cancies of different sizes, which act as buffers in the case of 
distortion of the skeleton as a consequence of its reaction 
with the non-skeletal cations. Considering the wide iso¬ 
morphism (mutual substitution) of cations in all crystallo¬ 
graphic positions (M, L, and X), the possibility of ions of dif¬ 
ferent degrees of oxidization and sizes occupying the skeletal 
positions, and also diverse variants of the skeletal charge 
compensation, the set of compounds and solid solutions, ei¬ 
ther known or predicted based on crystallochemical data, be¬ 
comes quite extensive. A controlled modification of the NZP 
phase compositions provides for ample diversity, smooth 
variations, and control of its properties. 

The increasing interest in materials with the NZP-type 
structure is due to the expanding possibilities of using them 
as high-tech ceramics that have the unique ability of not ex¬ 
panding under heating. This property of phosphates along 
with their high melting temperatures (Table 1) and physico- 
mechanical parameters (bending strength 1100 MPa, Young’s 


Fig. 1 . Crystal structure of NZP 
ceramics. The skeleton consist¬ 
ing of L-octahedra and X-tetra- 
hedra is indicated, as well as the 
nonskeleton positions Ml and 
M2. 
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TABLE 1 


Compound 

Melting 
(decomposi¬ 
tion) tempe¬ 
rature, K 

Compound 

Melting 
(decomposi¬ 
tion) tempe¬ 
rature, K 

LiZr 7 (P0 4 ) 3 

1963 

B a 0.5Zr2(PO 4 ) 3 

1973 

NaZr 7 (P0 4 ) 3 

1923 

B a 0.5Ti2(PO 4 ) 3 

1641 

KZr 2 (P0 4 ) 3 

1883 

Bai + v Zr 4 P 6 _ 2,Si 2x 0 24 

1773 

Na 5 Zr(P0 4 ) 3 

1433 

Na 2 FcZr(P0 4 ) 3 

1273 

Mgo.sZMPGtb 

1273 

L n o.33Zr 2 (P0 4 ) 3 

(Ln = La, Gd, Eu) 

1173 

C a o.5Zr2(P0 4 ) 3 

1973 

Zr 3 (P0 4 ) 4 

1173 

S r o.5Zr 2 (P0 4 ) 3 

1973 

M 0 3 3Nbi 67 (P0 4 ) 3 
(M = Mg, Mn, Co, Ni) 

1173 


modulus 50 GPa, Rockwell hardness HRA 90, and density 
3200 - 3650 kg/m 3 ), chemical and radiation resistance, and 
low thermal conductivity makes them suitable for applica¬ 
tion in products that should have high resistance to thermal 
shock (refractory lining and its elements, fittings for high- 
precision soldering, semiconductor substrates, catalyst carri¬ 
ers, optical benches), in electronics (ceramic electrolytes, 
bases), in machine building and transport subsystems (engine 
components), and in power engineering (radiation-resistant 
heat-insulating materials, localizing matrices for safe storage 
of radioactive waste). Note that one of the common names 
used for the variable-composition phases Na, + v Zr 7 P 3 _ v Si v O p 
with a NZP-like structure is NASICON. Other compounds 
having high ionic conductivity were named LISIKON, 
TITSICON, etc. Researchers of solid electrolytes who use 
these abbreviations relate samples to their composition and 
not to their structure. On the other hand, specialists in the 
sphere of ceramics dealing with materials that are structural 
analogs of the phosphate NaZr 7 (P0 4 ) 3 prefer the name NZP, 
which describes the family of compounds with the same type 
of structure and skeleton, although the overwhelming major¬ 
ity of this family are not sodium-bearing compounds and not 
ionic conductors. 

In view of the start of production of NZP ceramics [1] 
and expansion of its application areas, it becomes necessary 
to determining its general regularities relating its composi¬ 
tion and properties. The purpose of the present study is to de¬ 
termine the regularities correlating the composition, struc¬ 
tural transformations, and thermophysical properties of NZP 
materials, as well as the thermal resistance of ceramics and 
admissible fluctuations of temperatures and its gradients. 
These properties, such as thermal expansion, heat capacity, 
and thermal conductivity, have fundamental significance re¬ 
gardless of the application areas of these ceramics. 

The most characteristic property depending on the distri¬ 
bution of the force of the chemical bonds in a structure is 
thermal expansion: the less strong the bonds between the 
atoms, the greater the thermal expansion. For NZP-like lat¬ 
tices the thermal expansion is determined by the TCLE mea¬ 
sured parallel to a c and perpendicular to a a , which is the 


main crystal axis, and the average TCLE is a av = 
(2a fl + a c )/3. 

The majority of little-expanding materials of the NZP 
family [3, 4] have anisotropic thermal expansion and large 
absolute TCLE values and, nevertheless, low average TCLEs, 
for instance NaD 2 (P0 4 ) 3 (D = Zr, Hf, Ti), RbZr 7 (P0 4 ) 3 , 
Ca 0 5 Zr 7 (P0 4 ) 3 , Cd 0 5 Zr 2 (P0 4 ) 3 (Table 2).This is due to the 
fact that the TCLE of these particular compounds along the 
directions determined by the crystal symmetry have opposite 
signs. Similarly, the best of the known little-expanding mate¬ 
rials, such as zircon, cordierite, or quartz glass, have virtually 
zero expansion (a av = (0.5 - 4.2) x 10 6 K 1 in the tempe¬ 
rature interval of 273 - 1273 K). The consequence of differ¬ 
ent thermal expansion of different components in the volume 
of the article or in different directions in a single-phase mate¬ 
rials is mechanical stress on the grain boundaries, sufficient 
for causing cracking and disturbing the continuity of articles. 

At the same time, some NZP compounds have low 
thermal anisotropy and low TCLE values (less than 
2 x 10 ~ 6 K. 1 ). Such properties are found in phosphates 
CsZr 2 (P0 4 ) 3 , CsHf 2 (P0 4 ) 3 , CSj 3 Ln 0 3 Zr t 7 (P0 4 ) 3 

(Ln = Pr, Sm, Gd), Sr 0 5 Zr 2 (P0 4 ) 3 , Ca 0 5 _ v Sr v Zr 2 (P0 4 ) 3 , 
Sr 0 25 Zrj 5 Nb 0 5 (P0 4 ) 3 (Table 2). These phosphates, due to 
their ultralow thermal expansion, including near-zero TCLEs 
and anisotropic expansion, are capable of withstanding mul¬ 
tiple sharp fluctuations in thermal loads (thermal shocks) and 
are the most promising for developing thermomechanically 
stable materials. 

Based on our studies and the results obtained by other 
authors, we determined a system of empirical regularities of 
thermal expansion in NZP materials [4], which makes it pos¬ 
sibility to extrapolate these regularities to materials of the 
same structural type that have not yet been studied and per¬ 
form a justified selection of ceramic compositions with opti¬ 
mum thermal expansion parameters. 

When a NZP compound, for instance, NaZr 7 (P0 4 ) 3 is 
heated, the weakest bonds Na(Ml) - O expand much more 
than the strong bonds Zr(L) - O and P(X) - O; the crystal 
undergoes anisotropic expansion; it expands along the c axis 
and compresses along the a axis. As the result, the average 
TCLE of NZP compounds is very low. A series of com¬ 
pounds AD 7 (P0 4 ) 3 (D = Ti, Sn, Hf, Zr; A = Na, K) and 
LiGe 7 (P0 4 ) 3 exhibited a negligibly small variation of the 
size of the polyhedra D(L)0 6 and P0 4 with increasing tem¬ 
perature; therefore, the skeleton {[L 2 (P0 4 ) 3 ] - } 3oo of such 
compounds has thermal expansion approaching zero. Due to 
the arrangement of the less rigid (extendable) octahedra 
A(M1)0 6 in columns parallel to the c axis (Fig. 1), the heat¬ 
ing effect leads to the expansion of the structure in the direc¬ 
tion c (a c >0). In this case we observe a correlated rotation 
of the L0 6 octahedra and phosphor tetrahedra connected by 
their vertexes, and the structure compresses in the direction a 

K<o). 

A correlation of the effects of heating and crystallo- 
chemical substitution on the geometrical parameters of the 
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TABLE 2 


TCLE, l(r 6 K- 


Compound 


Aniso¬ 

tropy 

I a (i - aj 
l(r 6 K- 


Degree of 
anisotropy 

a„ -aJ 




Published 

source 


NaZr 2 (P0 4 ) 3 

-5.50 

22.30 

3.80 

27.80 

7.30 

[5] 

KZr 2 (P0 4 ) 3 

-5.30 

5.41 

- 1.73 

10.70 

6.20 

[5] 

RbZr 2 (P0 4 ) 3 

-8.00 

13.00 

- 1.00 

21.00 

21.00 

[6] 

CsZr 2 (P0 4 ) 3 

-0.56 

0.45 

-0.22 

1.00 

4.50 

[6] 

NaHf 2 (P0 4 ) 3 

-6.55 

21.00 

2.63 

27.50 

10.50 

[6] 

CsHf 2 (P0 4 ) 3 

- 1.10 

1.00 

-0.40 

2.10 

5.30 

[6] 

LiTi 2 (P0 4 ) 3 

0.49 

30.80 

10.60 

30.30 

2.90 

[7] 

NaTi 2 (P0 4 ) 3 

-5.30 

20.80 

3.40 

26.10 

7.70 

[7] 

KTi 2 (P0 4 ) 3 

-0.13 

6.62 

2.12 

6.70 

3.20 

[7] 

LiGe 2 (P0 4 ) 3 

-0.75 

22.50 

7.00 

23.20 

3.30 

[8] 

Ca 0 . 5 Zr 2 (PO 4 ) 3 

-2.57 

7.74 

0.87 

10.30 

11.80 

[9] 

S r o.5Zr 2 (P0 4 ) 3 

2.24 

2.28 

2.25 

0.00 

0.00 

[9] 

Cd 0 . 5 Zr 2 (PO 4 ) 3 

-3.50 

10.20 

1.10 

13.70 

12.50 

[5] 

Bao.5Zr 2 (P0 4 ) 3 

5.40 

- 1.80 

3.00 

7.20 

2.40 

[10] 

,375Zr 4 P5. 25 Si 0 75 0 24 

0.14 

0.46 

0.25 

0.32 

1.28 

[2] 

Sr 05 Hf 2 (PO 4 ) 3 

1.66 

1.67 

1.66 

0.00 

0.00 

[9] 

Cao.5Ti 2 (P0 4 ) 3 

6.99 

1.31 

5.10 

5.70 

1.10 

[9] 

S r o.5Ti2(P0 4 ) 3 

9.88 

-0.14 

6.54 

10.00 

1.50 

[9] 

Cao.25Sro.25Zr 2 (P0 4 ) 3 

-0.70 

1.10 

-0.10 

1.80 

18.00 

[H] 

S r o.25^ r i .5Nb 0 .5(PO 4 ) 3 

- 1.30 

2.60 

0.00 

3.90 

- 

[12] 

Na 3 CdZr(P0 4 ) 3 

7.04 

13.90 

9.30 

6.90 

0.70 

[5] 

, 3 Pr 0 3 Zrj 7 (P0 4 ) 3 

- 1.50 

0.50 

-0.80 

2.00 

2.50 

[6] 

Csj 3 Sm 0 3 Zr, 7 (P0 4 ) 3 

-0.90 

1.30 

-0.20 

2.10 

10.50 

[6] 

Csj 3 Gd 0 3 Zr j 7 (P0 4 ) 3 

-0.52 

-0.46 

-0.50 

0.10 

0.20 

[6] 

NbZr(P0 4 ) 3 

-3.40 

1.20 

-1.90 

4.60 

2.40 

[13] 

Mgo.33Nb! 67 (P0 4 ) 3 

-2.30 

3.40 

-0.39 

5.70 

14.60 

[14] 

Co 0 .33Nbj 67 (P0 4 ) 3 

- 1.50 

4.40 

0.47 

5.90 

12.60 

[14] 

Nio.33Nb] 6 7 (P0 4 ) 3 

-2.20 

3.00 

-0.47 

5.20 

11.10 

[14] 

Al 05 Nb L5 (PO 4 ) 3 

- 1.10 

~0 

-0.73 

1.10 

1.50 

[15] 

Fe 0 5 Nbi 5 (P0 4 ) 3 

-2.60 

1.30 

- 1.30 

3.90 

3.00 

[15] 


crystal lattice shows that with increasing tempe¬ 
rature an octahedron, for instance Na(Ml)0 6 , ex¬ 
pands but not enough for sodium to be replaced 
by a larger alkaline metal cation and obtain the 
polyhedron A(M1)0 6 that is equivalent in size 
(A = K, Rb, or Cs). Therefore, the crystallo- 
chemical substitution effect and the thermal ex¬ 
pansion effects have similar features: they are not 
equivalent, although they affect the same atomic 
bonds in the structure. 

Knowing the effect of the occupied and va¬ 
cant cation positions on the thermal expansion of 
a structure, one can predict the thermal expansion 
characteristics of NZP materials. The experimen¬ 
tal data accumulated suggest that Ml positions 
have a deciding effect on TCLE modification. 

In the case of compounds described by the 
general formula AD,(P0 4 ), the positions Ml are 
completely filled by cations A + and the positions 
M2 are vacant. If an occupied octahedron A0 6 is 
already extended along the crystallographic c axis 
at room temperature due to the introduction of a 
larger cation, the further expansion of the struc¬ 
ture along the c axis in heating will be suppressed 
by the bridge P0 4 tetrahedra that connect the 
skeletal framework columns. For instance, the ax¬ 
ial coefficients in the phosphates series 
AD 2 (P0 4 ) 3 (A = Na, K, Rb, Cs; D = Zr, Hf) in 
the temperature range from 293 to 1073 K are 
within the limits : a fl = (- 8 to - 0.56) x 10“ 6 ! , 

a c = (22.3 - 0.45) x 10~ 6 K _1 passing from the 
composition with Na (ionic radius r Na+ = 1.02 A) 
to the composition with Cs (r + = 1.67 A). 

TCLEs approaching zero and nearly zero aniso¬ 
tropy of the axial expansion are reached in the 
compounds CsD 7 (P0 4 ) 3 (D = Zr, Hf) (Table 2). Here a ge¬ 
neral tendency is observed: the larger cation filling the 
non-skeletal position Ml decreases both thermal expansion 
and the anisotropy of thermal expansion. The expansion - 
compression effects are the lowest in compounds with the 
smallest structure-forming cations D of the skeleton and 
grow as the size of the cation D increases. This relationship 
is determined by the fact that the more polarized metal ion to 
a greater extent localizes electron density on the bonds 
P - O - (D) and the double rotation of the P0 4 tetrahedron 
and D0 6 octahedron under a change of temperature becomes 
less facile. 

The presence of a 3D canal lattice in the NZP structure 
enables the migration of Li + and Na + cation from fully occu¬ 
pied Ml positions to vacant M2 positions as the temperature 
increases. This process is accompanied by more intense re¬ 
pulsion between the oxygen ion chains around the vacant Ml 
positions, which increases a c . Thus, for LiTi,(P0 4 ) 3 the dis¬ 
ordering of lithium cations with respect to positions Ml and 


M2 under heating leads to the maximum TCLE along the c 
axis ( a c = 30.8 x 10” 6 K 1 ) of all the NZP materials. 

If the filling of positions Ml is lower than unity and the 
positions M2 are vacant, which occurs, for instance, in com¬ 
pounds of the type B 0 5 D,(P0 4 ) 3 (B = Ca, Sr, Cd; Ba; D = Ti, 
Zr) or Ln 0 33 D 2 (P0 4 ) 3 , it is necessary to take into account the 
contribution of the cation-filled and vacant Ml positions to 
the relative expansion - compression of the structure with in¬ 
creasing temperature. It has been noted that the filling of Ml 
positions expands the structure along the c axis and com¬ 
presses it along the a axis. At the same time, vacant Ml posi¬ 
tions can either expand or compress in heating. In this case 
the electrostatic mutual attraction of the cation in the position 
Ml and the oxygen ion disappears, and the size of the vacant 
M1 position is controlled by the 0-0 contact depending on 
the value and the charge of the cations in the filled Ml posi¬ 
tion and in the skeletal framework. Thus, in the structures 
Cd 0 5 Zr 2 (P0 4 ) 3 and Ca 0 5 Zr 2 (P0 4 ) 3 the vacant positions Ml 
expand with increasing temperature and in Sr 0 5 Ti,(P0 4 ) 3 
they become compressed. By combining NZP compounds 
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TCLE, 10- 6 K- 1 



n 


TCLE, 10~ 6 1 



n 


Fig. 2. Dependence of TCLE a a (I) and a c (II) on the degree of 
filling M-positions in phosphates: a) Na 52i B v Zr(P0 4 ) 3 , where 
B = Mg (1 ), Ca (2), Sr (3); b) Na, _ NbZr,_ r (P0 4 ) 3 (4), 
Na , zr 2.25-0.25,( p ° 4 ) 3 (5), ^ l+ Zr 2 (V 3 _Si x )O l2 (6), 

Na, + Y Zr 2 _ x (P0 4 ) 3 (7), Na, +2 Mg Zr 2 _ x (P0 4 ) 3 (8). 


with TCLEs that are close in values and opposite in signs and 
can form wide ranges of solid solutions, it is possible to de¬ 
sign ceramic compositions with a controllable TCLE, includ¬ 
ing an extralow TCLE with minimal anisotropy and, accord¬ 
ingly, resistant to microcracking under sharp temperature 
fluctuations. For instance, ceramic materials based on the 
compounds Ca 05;t Sr t Zr 2 (PO 4 ) 3 and K 2v Sr 0 5 v Zr 2 (PO 4 ) 3 
have thermal anisotropy close to zero [11, 16]. This is due to 
the fact that KZr 0 ( P0 4 ), and Ca 0 5 Zr 2 (P0 4 ) 3 have positive 
TCLEs along the c axis (a c > 0) and negative values along 
the a axis ( a a < 0), whereas Sr 0 5 Zr 2 (P0 4 ) 3 has positive 
TCLE with (a a > a c ). Consequently, intermediate composi¬ 
tions with close to zero expansion anisotropy may exist be¬ 
tween the end members of the solid solution series. The re¬ 
sults of the research established that ceramic materials of 
the compositions K 0 2 Sr 0 4 Zr 2 (P0 4 ) 3 (a av = 3.1 x 10 _6 K _1 , 

I a fl - aj = 0.95 x 10 " 6 K~ 1 ) and Ca 025 Sr 025 Zr 2 (PO 4 ) 3 
(a av = - 0.1 x 10- 6 K->, | a a - aj = 1.8 x 10“ 6 k ') have 
minimal thermal anisotropy. 

If the positions M1 and M2 are not filled (with the skele¬ 
tal charge p = 0), the rotations of the polyhedra are less in¬ 
tense and the intensity of the thermal expansion of the lattice 
of the “empty” skeletal changes as the ionic radius of the 
atom in the L position changes. Under thermal loading the 
expansion - compression of the lattice along the c axis due to 
the presence of vacancies in the direction c depends on the 
strength of the 0-0 and L - O bonds. The average TCLE 
a av of the empty skeletons changes from - 1 x 10 “ 6 to 
1 x 10- 6 K- 1 . 

The crystallochemical nature of thermal deformations is 
more complex when the positions Ml are vacant and M2 are 
filled with cations. Thus, most compounds Na 3 R 2 (P0 4 ) 3 
(R = Sc, Ti, V, Cr, Fe, Y, In, Yb) undergo reversible phase 
transformations consisting in the adaptation of the structure 
to the changing exterior conditions (temperature, or specifics 


of the synthesis of samples). The existence of polymorphism 
in the compounds A 3 R 2 (P0 4 ) 3 (A = Li, Na) restricts the pos¬ 
sibility of their application as high-strength and crack-resis¬ 
tant materials. The data on the thermal expansion of com¬ 
pounds with totally filled M2 positions and vacant Ml posi¬ 
tions are limited to a few sodium-bearing phosphates. The 
ceramics based on these phosphates are highly-expanding 
materials (a av > (6 - 8) x 10“ 6 K~ *), since the average 
TCLE value of the oxygen polyhedra containing Na is high. 

Substitution with filling space Zr 4+ —» Na + + R 3+ in the 
crystal structure NaZr 2 (P0 4 ) 3 leads to the formation of solid 
solutions Na 1 + i R v Zr 2 _ x (P0 4 ) 3 , in which cations partly fill 
Ml and M2 positions. Three effects are observed in such 
substitutions with filling space: an increasing filling of M2 
positions, which causes the skeleton to expand along the a 
axis and compressing along the c axis; a change in the num¬ 
ber of vacancies in Ml positions leading to the expansion - 
compression of the Ml positions (vacancies); a very insignif¬ 
icant skeleton expansion due to the replacement of some Zr 4+ 
cations by larger cations. The dependences of the TCLE ( a a 
and a ( ,) on the degree of filling of the M-positions in the 
compounds and solid solutions investigated by us [5, 17] or 
described in the literature [13]: Na 5 _ 4v Zr T Zr(P0 4 ) 3 , 
Na 5 _ 2 BZr(P0 4 ) 3 (B = Mg, Ca, Sr), Na,_ 3 /eZr(P0 4 ) 3 , 
Na 3 _ v Nb A Zr 2 _ x (P0 4 ) 3 , Na 3 + A Zr 2 (P 3 _ x Si x )0 12 , 

Na l +.r Y v Zr 2-.v( P ° 4 ) 3 , and Na l+2z M gx Zr 2-*( P ° 4 ) 3 are COm - 

plex and nonmonotonic. 

With increasing degree of filling of the M-positions, the 
TCLE along the a axis (o . a ) is usually negative (Fig. 2); 
sometimes it grows and transforms from a negative into a 
positive one, whereas a c remains large (a c > 0). This ten¬ 
dency persists up to the filling M-positions equal to 1-2.5. 
The axial anisotropy in this case is the lowest. As the filling 
of M-positions grows from 2 - 2.5 to 3.5, the TCLEs for both 
directions may be positive, whereas a a usually decreases and 
the anisotropy increases again. This change in the TCLE may 
be explained as follows. Under a certain degree of filling of 
the positions Ml and M2, the thermal expansion effect along 
the a axis becomes comparable to the expansion effect along 
the c axis, and later on expansion along the c axis may be¬ 
come the predominant one. 

With full occupancy of Ml and M2 positions by cations, 
for instance in the phosphate Na,Zr(P0 4 ) 3 , the deciding ef¬ 
fect in thermal expansion parameters is exerted by the size of 
the cation in the M-positions and the average thermal expan¬ 
sion of the metal (M + ) - oxygen bond. One can expect a de¬ 
crease in the absolute values of a a and a c and a decrease in 
anisotropy as the size of the cations filling all M-positions 
grows. 

The above regularities of thermal expansion of NZP-type 
compounds provide a scientific interpretation of the pro¬ 
cesses of formation of materials with required properties, 
i.e., makes it possible to deliberately design new ceramic ma¬ 
terials with controlled TCLE, including a low TCLE with 
near-zero anisotropy. 
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The precision determination of the heat capacity that 
characterizes the quantity of thermal energy required for a 
certain change in the temperature of the material and thermal 
conductivity that characterizes the quantity of heat transmit¬ 
ted through the material under a preset temperature gradient 
are labor-consuming operations. That is why data on heat ca¬ 
pacity and thermal conductivity are limited to a few repre¬ 
sentatives of the NZP family. 

The heat capacity of different materials in the NZP fam¬ 
ily is 0.5 - 1.0 J/(g • K) in the temperature interval of 
273 - 623 K (Fig. 3). An important property of specific heat 
is that it is independence of the microstructure of the material 
(apparent density, porosity, distribution of pores by shape 
and size) at high temperatures. The heat capacity of most ma¬ 
terials investigated grows monotonically with increasing 
temperature and change insignificantly at temperatures 
above 500 K. The heat capacity curves of Na 5 Zr(P0 4 ) 3 , 
Na 3 Fe 2 (P0 4 ) 3 , and Na 3 Cr 2 (P0 4 ) 3 exhibit anomalies (phase 
transformations) related either to a modification of the posi¬ 
tion ordering of sodium cations in the skeletal vacancies of 
the first compound (transformation temperature 406.9 K) or 
to the processes of antiferromagnetic ordering for the two 
latter compounds (at temperatures of 47 and 12 K, respec¬ 
tively). 

Data on heat capacity and other thermodynamic func¬ 
tions are needed to study the thermodynamic aspect of syn¬ 
thesis and stability of NZP materials. The use of thermody¬ 
namic methods makes it possible to learn which of the possi¬ 
ble chemical reactions may occur spontaneously at the preset 
conditions (temperature, pressure, concentration) and how 
the conditions should be modified for the specified process 
to proceed in a required direction with a sufficient product 
yield. 

The initial materials in solid-phase synthesis of NZP ce¬ 
ramics are usually alkali metal salts, zirconium oxide, and 
ammonium dihydrophosphate (the chemical formulas of the 
compounds are followed by a letter in brackets indicating 
their physical state: “c” for crystal and “g” for gas): 

NaCl(c) + 2ZrOi(c) + 3NH 4 H 2 P0 4 (c) -> 

NaZr 2 (P0 4 ) 3 (c) + 3NH 3 (g) + 4H 2 0(g) + HCl(g); (1) 


Cp , J/(g ■ K) 



Fig. 3. Temperature dependences of heat capacity C p of crystalline 
compounds: Na 5 Zr(P0 4 ) 3 (1 ), Na 3 Cr,(P0 4 ) 3 (2), Na,MgZr(P0 4 ) 3 
(2), Na 3 Fe 2 (P0 4 ) 3 (4), Na 4 Zr 2 Si 3 0 12 (5), Na 3 Zr,Si 2 P0 12 ( 6 ), 
NaZr 2 (P0 4 ) 3 (7) and CsZr 2 (P0 4 ) 3 (8). 


tive reactants. The entropy (the measures of irreversible en¬ 
ergy dissipation) of the reactions A ( .5° (298) were calculated 
based on the values of the absolute entropy of the reactants at 
the temperature specified. The standard thermodynamic 
functions A r H°(T) and A r S°(T) at temperatures T above 
298.15 K were calculated based on general thermodynamic 
relations provided that the algebraic sum of the heat capaci¬ 
ties of the reactants taking into account their stoichiometric 
coefficients A (n C G p ) at 298.15 K is constant in the tempera¬ 
ture interval from 298.15 K to T: 

A,.H° (T) = A,.H° (298) + A (n C° p )(T - 298.15); 

A r S°(T) = A, 5'° (298) + A (n C° ) In (T/298.15). 

The standard Gibbs function of the reactions at all tem¬ 
peratures was calculated from the equation 

A ,.G°(T) = A r H° ( T ) - TA,.S° (T). 

If A,.G° <0, the reaction tends to be spontaneous in the 
direction “reactants —» products.” 


CsCl(c) + 2Zr0 2 (c) + 3NH 4 H 2 P0 4 (c) -A 
CsZr 2 (P0 4 ) 3 (c) + 3NH 3 (g) + 4H 2 0(g) + HCl(g); (2) 

5NaCl(c) + Zr0 2 (c) + 3NH 4 H 2 P0 4 (c) -A 
Na 5 Zr(P0 4 ) 3 (c) + 3NH 3 (g) + 2H 2 0(g) + 5HCl(g). (3) 

To analyze the conditions of the formation of 
NZP products, standard thermodynamic functions 
for their reactions of synthesis have been calculated 
(Table 3). The enthalpy (thermal effects) of the re¬ 
actions considered A / ./7°(298) at a temperature of 
298.15 K were calculated in accordance with the 
Hess law based on the standard formation enthalpy 
at the same temperature A ,H° (298) of the respec- 


TABLE 3 


Compound 

Cp (298), 

J/(g ■ K) 

- A r ZZ° (298), 

kJ/g 

-A,A 0 (298), 

J/(g ' K) 

kJ/g 

Published 

source 

NaZr 2 (P0 4 ) 3 

0.6238 

10.67 

2.25 

10.00 

[18] 

KZr 2 (P0 4 ) 3 

- 

10.43 

- 

- 

[19] 

RbZr 2 (P0 4 ) 3 

- 

9.59 

- 

- 

[19] 

CsZr 2 (P0 4 ) 3 

0.5122 

8.83 

1.88 

8.27 

[19] 

Na 5 Zr(P0 4 ) 3 

0.7919 

11.37 

2.36 

10.67 

[20] 

Na 3 Fe 2 (P0 4 ) 3 

0.7586 

- 

2.51 

- 

[21] 

Na 3 Cr 2 (P0 4 ) 3 

0.7570 

- 

2.62 

- 

[21] 

Na 3 MgZr(P0 4 ) 3 

0.7473 

- 

2.55 

- 

[21] 

Na 3 Zr 2 Si 2 POn 

0.6623 

10.96 

2.09 

10.34 

[22] 

Na 4 Zr 2 Si 3 0 12 

0.6632 

11.41 

2.08 

10.79 

[22] 
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X, W/(m ■ K) 



Fig. 4. Temperature dependences of thermal conductivity X of 
high-density (1-5) and low-density (6-9) types ofNZPceramics: 
C a o g Mgo iZ r 4 (P0 4 )6 (1 ), Ca Q9 Zr 405 (PO 4 ) 6 (2), Ca Q9 Li Q2 Zr 4 (P0 4 ) 6 
(3), Ca 08 Mg 02 Zr 4 (PO 4 ) 6 (4), Ca fls Mg 05 Zr 4 (P0 4 ) 6 (5), Zr 3 (P0 4 ) 4 
(6), CsZr 2 (P0 4 ) 3 (7), NaZr 2 (P0 4 ) 3 (S), and Na 5 Zr(P0 4 ) 3 (9). 

It can be seen from the data in Table 4 that the reactions 
used for the synthesis of NZP materials are endothermic. The 
Gibbs functions at 298.15 K of the reactions of synthesis are 
positive, i.e., the equilibrium at the temperatures specified is 
shifted toward the initial materials. The synthesis tempera¬ 
tures above which the processes (1) - (3) are thermodynami¬ 
cally resolved at a standard pressure are not high and amount 
to approximately 430 K for NaZr 2 (P0 4 ) 3 , about 410 K for 
CsZr 2 (P0 4 ), and about 570 K for Na 5 Zr(P0 4 ) 3 . For this rea¬ 
son, ceramic technology is a common method for producing 
NZP materials. 

To predict application areas for NZP ceramics and im¬ 
prove their production technologies, one has to know the rate 
of temperature change in the article under heat treatment or 
during service, i.e., its thermal conductivity and heat capac¬ 
ity. The thermal conductivity of ceramics depends both on its 
crystal structure (the more complicated is the crystal lattice 
and the more defective its structure, the lower is the thermal 
conductivity) and its ceramic structure (organization of crys¬ 
tallites of the same or different compositions in the volume 
of polycrystalline bodies). The temperature dependences of 
the thermal conductivity of known high-density (porosity be¬ 
low 10%) and low-density (porosity 26 - 38%) kinds of NZP 
ceramics are shown in Fig. 4. It can be seen that the 
dependences within the temperature interval considered have 
no peculiarities and the thermal conductivity grows smoothly 


with increasing temperature. The dependence of the thermal 
conductivity of NZP ceramics on its porosity is obvious in 
the physical context. As the thermal conductivity of low-den¬ 
sity ceramics is an intermediate value between the thermal 
conductivity of the main material and the thermal conductiv¬ 
ity of air contained in the pores between the particles, it is 
evident that due to the large difference between the thermal 
conductivity of the solid and gaseous phases, the looser ce¬ 
ramic with a higher content of air and with intense distur¬ 
bance of intergrain contacts will have the lower thermal con¬ 
ductivity. For low-density materials the thermal conductivity 
with zero porosity can be calculated from the formula [23] 

x=x l + 2E p (l-Z 0 /Z a )(2Z 0 /Z a +l)- 1 
° l-F p (l-Z 0 /Z a )(Z 0 /Z a + l)- 1 

where X, X 0 , and A a are the thermal conductivities of the po¬ 
rous body, its solid phase, and the gas (air) in its pores, re¬ 
spectively; V p is the volume fraction of pores in the porous 
body. 

The calculated thermal conductivity values with zero po¬ 
rosity within the temperature interval of 298 - 673 K are 
equal to (W/(m • K)): 0.6- 1.1 for Zr 3 (P0 4 ) 3 , 0.6-0.9 for 
NaZr 2 (P0 4 ) 3 , 0.7-1.1 for CsZr 2 (P0 4 ) 3 and 0.4 - 0.7 for 
Na 5 Zr(P0 4 ) 3 [24]. These values satisfactorily agree with the 
data on the thermal conductivity of high-density NZP cera¬ 
mics (A. = 0.6 - 1.4 W/(m • K), 298 < T< 873 K) [25] (U.S. 
patent No. 5102836). Thermal conductivity in thermal calcu¬ 
lations is usually taken as the mean working temperature of a 
particular product (material). 

It is notable that despite the spread in thermal conducti¬ 
vity values found in the literature caused by the chemical and 
phase compositions of the tested NZP materials, specifics of 
their crystal structure (valence of ions, their distribution in 
the crystal lattice, type and number of defects, the type of 
chemical bond), microstructural organization (size, disper¬ 
sion, shape of grains, and porosity type), as well as specifics 
of measuring methods, NZP ceramics has lower thermal con¬ 
ductivity than industrial stabilized zirconium dioxide 
(1.95 - 2.44 W/(m • K)) in the temperature interval of 373 - 
1673 K), which is the main component in widely used 
refractories for heat-enclosing structures of high-temperature 
furnaces. Ultralow thermal conductivity, high melting tem¬ 
peratures, ultralow thermal expansion with nearly zero aniso- 


TABLE 4 


r,K 

A r H°, kJ/g, for reaction 

\S°, J/(g 

• K), for reaction 

A r G°, 

kJ/g, for reaction 

(i) 

(2) 

(3) 

(1) 

(2) 

(3) 

(1) 

(2) 

(3) 

298.15 

1.061 

0.802 

1.679 

2.484 

1.990 

2.964 

0.320 

0.208 

0.796 

400 

1.058 

0.799 

1.674 

2.476 

1.983 

2.947 

0.068 

0.006 

0.495 

450 

1.057 

0.798 

1.671 

2.474 

1.980 

2.941 

- 0.007 

- 0.093 

0.347 

500 

1.055 

0.797 

1.668 

2.469 

1.978 

2.935 

-0.180 

-0.192 

0.201 

550 

1.054 

0.796 

1.665 

2.467 

1.976 

2.930 

-0.303 

-0.291 

0.054 

600 

1.052 

0.794 

1.663 

2.464 

1.974 

2.925 

- 0.426 

-0.390 

- 0.092 
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tropy, and chemical and radiation resistance suggest that 
NZP ceramics are promising heat-insulating materials ex¬ 
ceeding the known high-temperature insulators [23] in their 
capacity to withstand thermal shocks. 

Considering the growing interest in the industrial pro¬ 
duction of NZP ceramics that have a unique set of physical 
and physicochemical properties fortunately combined in the 
same material and using the results of analysis of its thermo¬ 
physical properties, taking into account crystallochemical 
data and chemical compositions will make it possible to de¬ 
velop new types of ceramics capable of resisting thermal 
shocks and possessing the required level of thermophysical 
parameters. 

The study was performed with financial support of the 
Russian Fund of Fundamental Research (project 02-03-32181). 
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